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Abstract 

Moduli, modulini and the gravitino have gravitational-strength interactions, 
and thermal collisions after reheating create all of them with roughly the same 
'. abundance. With their mass of order 100 GeV, corresponding to gravity-mediated 

I supersymmetry breaking, this leads to the well-known bound 7TR < 10^ GeV on 

■ the reheat temperature, where 7 < 1 is the entropy dilution factor. The vacuum 

O . fluctuation also creates these particles, with abundance determined by the solution 

of the equation for the mode function. Taking the equation in each case to be the 
one corresponding to a free field, we consider carefully the behaviour of the effective 
mass during the crucial era after inflation. It may have a rapid oscillation, which 
does not however affect the particle abundance. Existing estimates are confirmed; 
the abundance of modulini and (probably) of moduli created from the vacuum is 
• less than from thermal collisions, but the abundance of gravitinos may be much 

bigger, leading to a tighter bound on Tr if supersymmetry breaking is gravity- 
mediated. 
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Introduction Long-lived particle species constrain models of the early Universe, be- 
cause they are cosmologically dangerous. In this note we are concerned with species that 
have gravitational-strength interactions. According to present ideas, these include the 
gravitino (spin 3/2), and species usually called moduU (spin 0) and modulini (spin 1/2). 

The abundance of a given species is conveniently specified by the ratio n/s, where n is 
the number density and s is the entropy density, evaluated at the epoch of nucleosynthesis. 
A particle with a gravitational-strength decay rate, and mass 100 MeV to 10^ GeV, decays 
before the present but after nucleosynthesis. To avoid interfering with nucleosynthesis 
its abundance must satisfy ^ 

-<io-''. (1) 

s 

(The precise bound lies between 10~^^ and 10~^^, depending on the mass.) If supersym- 
metry breaking is gravity-mediated, the masses of the gravitino, moduli and modulini 
are all expected to lie in this range. 
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A particle with gravitational-strength interactions and mass < 100 MeV is stable 
(lifetime longer than the age of the Universe). This is expected to be the case for moduli, 
modulini and the gravitino, if supersymmetry breaking is gauge-mediated. Any stable 
species which is non-relativistic at the present epoch has matter density 

l]x = 4.4xl0^/.-2-^!^<0.3, (2) 
1 GeV s 

where h ~ 0.65 is the Hubble constant in units of lOOkms"^ Mpc~^. 

Gravitinos are created by particle collisions after reheating, with abundance at nu- 
cleosynthesis 

10-^^ (3) 



thermal lO^ GcV 

In this expression, Tr is the reheat temperature, and 7^^ > 1 is the increase in entropy, 
per comoving volume, between reheating and nucleosynthesis. Moduli and modulini are 
created with roughly the same abundance. If Eq. (|l]) holds for moduli, modulini or the 
gravitino, we arrive at the famous bound 7TR < 10^ GeV. 

Creation from the vacuum The vacuum fluctuation during inflation provides an- 
other mechanism for creating particles. Indeed, although there are no particles during 
inflation, each field has a vacuum fluctuation. During some era containing the end of 
inflation, the evolution of the mode functions may become non-adiabatic so that particles 
are created 0. After creation stops, the number density is 

n = n^{a^/af', (4) 

where a is the scale factor of the Universe, and a* is its value at the end of inflation. We 
shall refer to the constant as the number density at the end of inflation, even though 
particle creation stops only some time later. 

If the particles have only gravitational-strength interactions, or are very heavy 
they will not thermalize at reheating. The particles created by the vacuum fluctuation 
retain their identity, and ||^ their abundance at nucleosynthesis is 



n 
s 




(5) 

To^j • 

In this expression, H^, is the practically constant value of the Hubble parameter during 
inflation (to be precise, the value at the end of inflation), and V = 3MpH^ is the 

inflationary potential. The observed cosmic microwave background anisotropy requires 

1 

H if * < lO^'' GeV (K* < 10^^ GeV). The corresponding number density, for a stable 
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non-relativistic species mass mx, is 



- 10^°|^l(7^)l^|7 (7) 

\f (8) 




1014 GeV/ VlOi^GeV/ VlO^GeV 



If < lOiJ^, the thermal mechanism definitely dominates [0, but if is bigger the 
vacuum fluctuation may dominate, leading to a tighter constraint on Tr. In terms of the 
occupation number \l3k\^ for states with momentum fc/a, the number density per helicity 
state of a given species is 

ivr^ a-^ Jo 

Significant particle creation will occur only up to some fcmax, such that \Pk\ ~ 1 at 
k = fcmax, and < 1 at bigger k. 

For fermions, < 1 leading to the rather firm estimate P| 

~ 10-2(A;^,^/a,)3. (10) 

For scalar particles, this expression gives in general only a lower bound on n*, because 
\Pk\ can be much bigger than 1, but as we shall discuss it is likely to give very roughly 
the right order of magnitude for moduli. 

For a given particle species, fcmax is determined by the relevant mode function equa- 
tion. If we are dealing with a free field, this equation corresponds to the Klein-Gordon 
(spin 0) Dirac (spin 1/2) or Rarita-Schwinger (spin 3/2) equation, plus constraints to 
eliminate unwanted degrees of freedom except in the case of scalars. The resulting equa- 
tions for spin zero and spin 1/2 are well known, and they have recently been given also 
for helicity 3/2 gravitinos |], ^ and for helicity 1/2 gravitinos g, ^. 

In the following we consider each case in turn, paying special attention to the be- 
haviour of the effective mass that appears in the mode function equation. It may oscillate 
rapidly just after inflation, but this does not affect the estimate of fcmax- In agreement 
with existing estimates, we flnd that A;max ~ a^fn^:, where is the typical magnitude of 
the effective mass soon after the end of inflation. Also in agreement with these estimates 
we flnd that for moduli and modulini, m=K ~ H^, while for helicity 1/2 gravitinos ~ M 
where M > is the mass of the fleld responsible for the inflaton potential.^ 

^Throughout this paper, we assume that oscillation of the field responsible for the inflaton potential 
is dominated by its the mass term. This excludes the case of inflation with a quartic potential. In that 
case, the abundance of particles created by the quantum fluctuation is enhanced, being the same as it 
would be in the case of instant reheating jsj. In the case of the gravitino, we find that the advertised 
estimate is reduced in the unlikely event that the gravitino mass (appearing in the Rarita-Schwinger 
equation) is much bigger than H around the end of inflation. 
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Scalar particles The case of scalar particles, in particular moduli, is more complicated 
than the case of spin 1/2 particles. The basic reason is that the occupation number 
can be much bigger than I.Asa result the number density Eq. (^) may be dominated by 
long wavelength modes. Related to this is the fact that the classical scalar field can have 
a practically homogeneous value, which in the early Universe is different from the vacuum 
value. In that case, the field eventually oscillates about the vacuum value, corresponding 
to the existence of particles with negligible momentum. 

For a given field 0, a crucial quantity is the effective mass-squared m^{t) in the early 
Universe. (We take the origin to be a fixed point of the symmetries so that V = Vo+^m'^(j)'^ 
plus higher terms.) Assuming supergravity, it is typically of order ±iJ^, until H falls 
below the true mass mx- Consider first the value during inflation |]lT], |12|, 0. Unless it 
is dominated by the D term, the potential V is of the form exp(i^/Mp)X, where K is 
the Kahler potential, and the form of X does not concern us for the moment. We are 
assuming canonical normalisation, which means that K = to leading order. As a 
result the potential gives a contribution of the form 

ml, = V/M^ + ... . (11) 



The other contribution depends on X, and is of order ±V"/M|. In general there is no 

^pot 



reason to expect a cancellation, so we expect m'^^^ = XV /Mp, with A ~ ±1. This is the 



only contribution to during inflation, giving a practically constant value m'^^^ = 3XH^. 

Now consider the era just after inflation. The fleld x responsible for V oscillates, with 
angular frequency M equal to its mass. As a result V and m'^^^ oscillate with the same 
frequency. In addition, there will be a contribution [|14|, |^ '^kin ~ ^'x^/^p with A' of 



order ±1 but unrelated to A. During this era, the energy density p = SH'^M^, and the 
pressure P are given by 

P = V + \x' (12) 

p = -y+lx' (13) 

The quantities V and oscillate out of phase, so that does not oscillate. In fact, 
energy conservation gives H / = — 1(1 + P/p), and according to fleld theory \P\ < p. 
If A and A' were equal, m?{t) would have the same behaviour, but as they are not it 
oscillates. Redeflning the couplings we have 

c^H'^{\ + X' cos(Mt)) , (14) 

with |A| ~ I A' I ~ 1. 

The energy in the oscillation of x "will sooner or later be spread among several particle 
species, through reheating or earlier during preheating (through parametric resonance). 
Adding together the incoherent equivalent fleld oscillations, we shall now get a non- 
oscillating contribution m?{t) ~ ±if^. 
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Assuming that ~ ±iJ^ the cosmology of a scalar field depends crucially on the 
sign. If it is negative, the classical field value will be driven away from the fixed point, 
and will end up oscillating about its vev, corresponding to particle production that has 
nothing to do with the vacuum fiuctuation. It is this case that is generally expected 
to occur for, at least some, moduh. In the case of gravity-mediated supersymmetry 
breaking, the estimated abundance of moduli is cosmologically dangerous. This is the 
usual version of the moduli problem If, on the other hand, is positive during 



and after infiation, the classical value is always at the fixed point of the symmetries, and 
particle creation from the vacuum fiuctuation becomes relevant. This case, corresponding 
to the field held at a fixed point of the symmetries, may occur for at least some of the 



moduli as discussed for instance in [|I6|. We are trying to see if the moduli problem 
occurs also in this case. 

Before continuing with the case ~ H^, let us note that during infiation, it may 
be that \m?\ <^ H^. (This is actually mandatory for the infiaton field |]T^, which is not 



however expected to be a modulus.) If that happens, then (assuming minimal coupling 
to gravity) the very long wavelength quantum fiuctuation acts as a classical field dis- 
placement, which again eventually oscillates giving a situation essentially the same as 
the one in which < -H"^ |, |19|.|] 

Now we briefiy recall the formalism for the vacuum fiuctuation (see for instance 0). 
A suitably chosen mode function obeys the equation 

u" + w^M = , (15) 

where 

= k'^ + {amf (16) 



a" 



rfi^ = - (1 - 60 -T , (17) 



a 



and ^ = for minimal coupling to gravity. In these expressions the prime denotes 
differentiation with respect to conformal time t], where dt] = dt/a. (An overdot will 
denote d/dt.) Einstein gravity gives 

^ = Ih\1 - 3P/p) . (18) 

During infiation, a" /a? = 2H^, and during matter domination a" /o? = \H^- 

The mode function takes on a simple form during any era in which uo satisfies the 
adiabatic condition. This is the condition that u be slowly varying on the conformal 
timescale u!~^, or to be precise the conditions 

|cj'|<cj2, |cJ"|<cJ^ (19) 



^This conclusion does not assume that m? /H^ is small also after inflation, which can hardly occur 
since each particle species present gives a separate contribution to it. The calculation of |^ reaching an 
apparently different conclusion did not take into account the very long wavelength fluctuation. 
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During any era in which the adiabatic condition is satisfied, approximate solutions of 
Eq. (|15D are m'"^''' and u^'^''^*, where 

^adia ^ ^-i exp(-i r todr]) , (20) 

and ?7o is arbitrary. (Over a time interval -C H^^, these solutions reduce to the flat 
spacetime solutions.) 

The adiabatic condition is satisfied at early times during infiation, and one adopts 
the initial condition u = u^'^^'^, corresponding to the vacuum. At late times the adiabatic 
condition is again satisfied and 

u = au'^"' + Pu'^'^*, (21) 
with |ap — = 1. The occupation number is and an explicit formula is 

4uj\f3\^ = \u'\^ + uj^\u\^ -2uj. (22) 

There is negligible particle creation if <^ 1. A brief failure of the adiabatic 
condition need not lead to significant particle creation. At least for practical purposes, 
the criterion for significant particle creation is the failure of 

PI < cj^ , (23) 

where the average is over a conformal time interval uj~^. We shall call this the weak adi- 
abatic condition. It says that the fractional change in u is small, in a typical conformal 
time interval uj~^. To be convinced of its relevance, one can imagine the numerical inte- 
gration of Eq. (p!5D . Alternatively, one can note that Eq. (plSj) describes a two-dimensional 
harmonic oscillator with a time-varying frequency; the initial condition corresponds to 
circular motion, and the final condition will also correspond to circular motion if the 
weak adiabatic condition is satisfied. 

We want to estimate fcmax; defined as the wavenumber above which particle creation is 
negligible. Equivalently, fcmax is the value of k, above which the weak adiabatic condition 
Eq. ( p3D is satisfied. We assume |^| < 1, ~ H"^, and ~ i?^. If we took m? / H'^ to 
be constant, there would be no need for the average in Eq. (|2^) . Since aH peaks at the 
end of infiation, and varies on the inverse timescale H, this would give /Cmax 

~ a^H^ as 

advertised. In fact, changes at the end of infiation on the inverse timescale M >^H. 
If we ignored the average, this would give /Cmax ~ a^^/H^M. Restoring it, we have 
uj' = {H^/M)uj', leading again to fcmax ~ ct*-f^*- 

We still have to check that Eq. (P) is dominated by the upper limit, giving ~ 
10~^H^. Three independent calculations I^O], ^ ^ verify that this result is roughly correct 
if fh? /H"^ is fairly close to 1 during infiation.^ On the other hand, n^, is much bigger if 
m^/H^ is significantly below 1 ||^, n], due to the long wavelength effect mentioned earlier. 



^Reference pOl took |m| ^ H^, and ^ ~ 1, and solved the mode equation assuming a cx exp(_ff*t) 
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The conclusion is that the abundance of moduh created from the vacuum is less than 
the abundance from thermal collisions, if m?/H'^ is fairly close to 1 during inflation, 
which (in magnitude) is the 'best guess' from supergravity with the potential dominated 
by the F term. On the other hand, the abundance of moduli created from the vacuum 
is enhanced if m? / H"^ is significantly below 1 during inflation. Related to this is the fact 
that a negative value m? < —H^ (during or after inflation) generates moduli through the 
oscillation of the classical field, with abundance typically much bigger than the abundance 
from thermal collisions. 

Before leaving the subject of scalar particle creation, we should mention the case of 
a stable particle with true mass mx ^ 10^^ GeV [0, ^. Such a particle is a dark matte 
candidate (Wimpzilla) which is never in thermal equilibrium, and whose abundance from 
thermal collisions may be strongly suppressed. If mx accidentally happens to be of order 
if*, one again has fcmax ~ a*H^, leading again to ~ 10~^if^ if Eq. (^ is dominated by 
the large k limit. Using a variety of approaches to the solution of the mode equation, it 
was found that ~ 10~^if^ (taking ,^ = |) in rough agreement with this estimate. 
For a suitable choice of one has (Eq. (^) a dark matter candidate. Alternatively the 
Wimpzillas can be created through (suppressed) thermal collisions, parametric resonance 
or from a first-order phase transition . 

We end this discussion of scalar particle production by mentioning the exceptional 
case ^ = 1/6 with m — > 0. The adiabatic condition is then satisfied at all times and there 
is no particle production, but there is no reason to expect this case to occur. The lack of 
particle production reflects the fact that the field equation is invariant under conformal 
transformations of the metric, making the expanding universe indistinguishable from fiat 
spacetime. 



Spin 1/2 particles The formalism for the creation of spin 1/2 particles is described 
for instance in |^. A suitably normalized two-component mode function satisfies 




—am k 
k am 




(24) 



The evolution is unitary so that one can choose 

+ |m„|2 = 1. (25) 

during inflation, and either a oc fs (matter domination) or a cx ^2 (radiation domination) afterwards, 
with u and u' continuous at the transition. Reference took m = ChH (with Ch a constant of order 
1) and ^ < 1, and solved the mode equation by the same matching procedure, finding n■^, ^ 1Q~^'^H^. 
Reference made the same assumptions about m and ^, but solved the mode equation numerically 
assuming a quadratic inflation potential. The flrst two calculations ignore the oscillation, after inflation, 
of m and P/p- The third includes the latter oscillation but still ignores the former. 
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Defining u = u+ one findsQ 

u" + (u^ + i{may) u = , (26) 

with 

u^ = k^ + {maf . (27) 
If am is slowly varying, solutions of the two-component mode equation are 



u± = M^'^'^ = N± exp (-i f ujdri 



u± = 
where 

1 

The ratio / comes from Eq. (p^ , and the normalization from Eq. (pSj) . Over a time 
interval ^ -f/^^^, these are the usual flat spacetime solutions. The adiabatic condition, 
ensuring that they are valid, is 

(ma)' < uj^ , [ma)" < uj^ . (30) 



It may be obtained by inserting the solution into Eq. ( pq) ||21| . Equivalently, one can 
note that Eq. (^^ requires \N^/N±\ -C uj, while the second differentiation to get Eq. ( ^6]) 
requires also \N'J_/N±\ <^ cu^; the conditions on are the strongest, and lead to Eq. (|30|). 

The adiabatic condition is satisfied at early times during inflation, and at late times 
after inflation. Taking u± = u'^^'^ at early times, corresponding to the vacuum, one finds 
at late times 

M± = auf^ + pvf^ , (31) 
with |ap + = 1. This corresponds to 



u 



aN+ exp (^-i ^ cudr]^ + (3N_ exp [i ^ udi]^ . (32) 

An explicit formula for the occupation number is 

|/5|2 =(^^^am- 2Im(MM*')) / (2cu) . (33) 

Analogously with the scalar brief failure of the adiabatic condition will not 

cause significant particle creation, the condition for that being 



|(ma)'|<cj% (34) 



■*In (ma)' is incorrectly replaced by ma' . With the ansatz m/ H =constant that was used there, 
one can check that this has no effect if inflation is followed by radiation domination, and is equivalent 
to the replacement m/H — > ^m/H if inflation is followed by radiation domination. 
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where the average is over a conformal time interval uj~^. 

The adiabatic condition is satisfied in the hmit m ^ 0, refiecting the invariance of 
the Dirac equation under conformal transformations of the metric. There is no particle 
creation in this limit. 

In contrast with the case for scalars, spin 1/2 particles with gauge interactions usually 
have effective mass |m| <^ H m the early Universe [|1^. There is essentially no creation 
of such particles from the vacuum. 

For gauge singlets, including modulini, the situation is similar to that which we dis- 
cussed already for scalars and spin 1/2 particles. Supergravity is expected to generate an 
effective mass m ~ ±i?, coming from several different sources just as for scalar masses 
|T^. The ratio m/H is constant during inflation, and oscillates afterwards with angular 



frequency M >^ H^, before settling down to a new constant value. Applying the weak 
adiabatic condition Eq. (0) (not the full condition Eq. (|30D , which would give too big 
a valueQ for fcmax) this leads to the advertised estimate fcmax ~ and n* ~ 10~^H^, 

leading to the conclusion that the creation of modulini from the vacuum is less efficient 
than the thermal creation. 

Note that, in contrast with the scalar case, the abundance of spin 1/2 particles is not 
altered ii m/H happens to be very small during infiation, which may happen in some 
models. 

The gravitino mode equations At the time of writing, the study of gravitino pro- 
duction is just beginning 0, ||, ||, [l^. One has to consider separately the production 
of helicity 1/2 and helicity 3/2 particles, as seen by a comoving observer. If the grav- 
itino in the early Universe can be treated as a free field, it satisfies the Rarita-Schwinger 
equation, with time- dependent effective massP] 



m 



M/^'^W\/Ml, (35) 



where W is the superpotential (a holomorphic function of the complex scalar fields 0„) 
and K is the Kahler potential (a real function of 0„ and 0*). In minimal supergravity, 
K = I] |0nP but one expects additional terms suppressed by powers of Mp. The mode 
equations are found by imposing constraints to eliminate unwanted degrees of freedom. 
For helicity 3/2, one finds |, |10| the same equation, Eq. (H), as in the spin 1/2 



case. 

The mode equation for helicity 1/2 gravitinos is more complicated [|, 0. A suitably 



^In a different context ||2l| , fcmax is estimated on the assumption that there is significant particle 
production whenever the full adiabatic condition Eq. (|30| ) is violated. We have checked that, in that 
case, the same result would have been obtained using the more appropriate weak adiabatic condition 
Eq. 11). 

^This is the well-known expression of = 1 supergravity. In the early Universe, it may be reasonable 
to consider instead E3l iV > 1. 
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defined two-component mode function satisfies 



—am kG 
kG* am 



The function G may be written 
wfiere A = Ai + iA2 is the quantity defined in 

A, = 



G = e-'^A* 



{P/p)H^ - m^ 



(36) 
(37) 



A, 



^2 ^2 

2 m 

3 + ^2 ' 



(3^ 



and 



cos^^ 
sin 9 



m 



m^ 



2mH 
m^ _l_ ff2 



(39) 



The evolution defined by Eq. (|36D is unitary, leading to Eq. (25). The function 



u = G 2M_,_, satisfies 

with 

and 



u 



lima] 



u 



[ma] 



m 



1 3 , 3 , I A 
— m m/li HAo 

2 2 2 2zy4 

1 G 
"^+2^G 

m + -(argG) + 2^ 



(40) 
(41) 

(42) 

(43) 
(44) 

]; they are equivalent 



The first expression is given in [§, while the second is given in 
by virtue of Eq. (g3)-[| 

^To be more precise, the last terms of Eqs. ( ^0| ) and ( ^2[ ) appear with opposite sign in The origin 
of this discrepancy, which may be partly due to different conventions, is not clear and will not affect 
order-of-magnitude estimates of the abundance. References ^ give Eqs. (^-(^ only for the special 
case \A\ — 1, which may be the only one of physical interest, but we have checked that they follow from 
Eq. (H) also in the case |A| 7^ 1. 
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The gravitino abundance and its cosmological consequences The formalism 
of the last section applies if the gravitino can be treated as a free field in the early 
Universe, corresponding to the Rarita-Schwinger equation. This has been demonstrated 
IP] only when the energy density and pressure are dominated by the oscillation of a single 
homogeneous scalar field, with minimal kinetic term and a potential (taken to be tree- 
level) dominated by the F term of = 1 supergravity. It has been shown [§, ^ that in 
this case, \A\ = 1, so that the production of helicity 1/2 gravitinos is exactly the same 
as production of a spin 1/2 particle with mass m. At least if the magnitude of the scalar 
field is small on the Planck scale, this particle can be identified with the goldstino of 
global supersymmetry p, ||. 

In the case \A\ = 1, the last term of Eq. (^21) may be written 



2iA 2A2 2 [1 -w2 + 2(l + w)mV/72]2 

where w = P/p. During slow-roll inflation, w is close to —1 with {wl -C H(l + w)^, and 
afterwards w will oscillate between ±1, with angular frequency equal to the mass M >^H^ 



of the field responsible for the inflaton potential. Bearing all this in mind, Eq. (45) shows 
that if m < if, m increases sharply just after inflation, to a value of order M, in a time 
of order . Since we are dealing with the spin 1/2 mode equation, the adiabaticity 
condition is (am)' <C a;^, leading to the estimate fcmax ~ and ~ 1Q~'^M^ . 

This is the result already anticipated [§, ||. It is indeed reasonable to suppose that 
m <H, because according to supergravity if^ = ip(t)/M| always receives a contribution 
—iv?{t). At the present epoch this contribution is cancelled to high accuracy for some 
unknown reason (part of the cosmological constant problem) but during inflation there is 
no reason to expect that to happen, in the usual case that ii* is bigger than the present 
value of the gravitino mass. 

This estimate of the gravitino abundance is based on Eq. (^OD for the mode function, 
taking \A\ = 1 so that the equation becomes the same as in the spin 1/2 case. The 
equation |y4| = 1 (with A defined by Eq. (|38D) cannot be exact, because it relates the 
gravitino mass m{t), given by Eq. (|5BD, to the density and pressure. After reheating, 
or significant preheating, the density and pressure have nothing to do with the scalar 
fields that define m. One can even construct a model of inflation, albeit an unrealistic 
one, in which m{t) can be chosen independently of the density and pressure throughout 
the history of the Universe; this is the model |2^ which generates the inflaton potential 
entirely from a D term. 

It may be that whenever \A\ is different from 1, Eqs. (^) and (|40|) cease to be true. 
Let us nevertheless discuss briefly the opposite possibility, that the gravitino in the early 
Universe is described by the Rarita-Schwinger equation (plus constraints) with \A\ 7^ 1. 
This equation has well-known problems in a generic spacetime but it is not clear 
that they exist in the special case of the Robertson- Walker Universe. The phase velocity 
\A\ associated with Eq. (^) will be less than 1, as is presumably required by causality, 
provided that m does not vary too rapidly. Coming to the question of gravitino creation. 
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the following argument suggests that the estimate fcmax ~ a*M is likely to survive. 
Assuming as always m < H, it is clear that \m\a will have the same general behaviour 
as in the \A\ = 1 case, being of order a*M just after inflation and smaller at earlier 
and later epochs. Also, \A\ is at least roughly of order 1 just after inflation, since P/p 
oscillates between ±1. Choosing k > a^M, the mode equation therefore reduces to 



u" + uj'^u = 0, with u = \A\k. This is the same as Eq. ([T5|), and if Eq. (p3D is satisfied the 
positive-frequency mode is never generated, and there is no particle creation. Eq. (|23|) is 
just 

k>jA\'/\^\''- (46) 



It follows that fcjnax is of order a*M or the maximum of whichever is bigger. 

By inspection of Eq. (|38|), one can see that the latter quantity is typically of order a^^M 
just after inflation and smaller at earlier and later times, so that indeed fcmax ~ a*M. 

Finally, let us look at the cosmological significance of the estimate ~ 10~^M'^. 
Depending on the model of inflation, M can be anywhere between H^, and (roughly) 
V^^^, where V^^^ ~ {MpH^)^ is the inflaton potential. In the extreme case M ~ V^, 
Eq. (D gives § 

-~ 10-^-47 • 47 
s Vi 

In the case of gravity-mediated supersymmetry breaking, n/s < 10^^^ leading to a tight 
bound on Tr, which becomes more stringent as the scale of inflation decreases, with the 
nucleosynthesis bound Tr > 10 MeV imposing the constraint 

> 10^ GeV . (48) 



Conclusion We have considered the abundance of moduli, modulini and gravitinos 
created from the vacuum fluctuation, paying special attention to the behaviour of the 
effective mass just after inflation. In agreement with previous estimates, we find that 
creation from the vacuum is likely to be insignificant for moduli and modulini, but may 
be very important for gravitinos. These estimates are based on the assumption that in 
the early Universe, the relevant equation is the one corresponding to a free field. 
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